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Abstract 

Based on a recent paper by Rothe and Schafer on compact binary systems [J. Math. Phys. 
51 , 082501 (2010)], explicit expressions for canonical center and relative coordinates in terms of 
standard canonical coordinates are derived for spinless objects up to second post-Newtonian ap¬ 
proximation of Einstein’s theory of gravity (the third post-Newtonan order expressions are available 
in form of supplementary data). The inverse relations, i.e. the dependence of the standard canon¬ 
ical coordinates on the canonical center and relative coordinates, are also given up to the second 
post-Newtonian approximation. The famous Pythagorean-theorem-type Lorentz-invariant relation 
between the system’s total energy or Hamiltonian squared, the rest energy or mass squared - solely 
depending on relative coordinates ~, and the total linear momentum squared are explicitly shown 
through second post-Newtonian approximation. 


* grgira001@myuct.ac.za 
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I. INTRODUCTION 


Compact binaries are important objects in the relativistic astronomy and astrophysics. 
Their orbital motions however are difficult to calculate if the Einstein theory of gravity is 
deeply needed for sufficient accuracy. Analytically, only approximate solutions are known 
whereby post-Newtonian (PN) approximation calculations have shown to be most efficient, 


e.g 


q: 


Very helpful is the knowledge of Hamiltonians because they easily allow the tran¬ 
sition to the rest frame by putting the total linear momentum, which simply is the linear 
sum of the linear momenta of the single components, equal to zero which is allowed when¬ 
ever the center-of-mass motion is unimportant. In this way, the rest-frame dynamics of 
post-Newtonian gravitationally bounded binary systems has often been derived, e.g. 2|-l6j. 
However, if for example a binary system is part of a hierarchical triple system its total motion 
becomes crucial and a more general treatment is required, e.g. ItHI^. The same holds if the 
recoil of merging binary systems through gravitational wave emission cannot be neglected, 
see e.g. jffi, 12|. A specihc topic in the gravitational recoil problem is the generalization 
of the effective-one-body method to the binary-center motion, see e.g. j^|. These topics 
clearly make a strong case for the investigations performed in the present paper. 

The transition from standard canonical coordinates to center and relative coordinates is 
nothing but a canonical transformation. The generators of this canonical transformation 
have been derived up to third post-Newtonian (3PN) order, even including spin degrees 
of freedom, in Q. Howeve. explic.t exp..io. tox the xe.atxons between the Centex and 
relative coordinates and the standard canonical coordinates of the binary components, and 
vice versa, have not been given. The present paper closes this gap through 2PN order for 
spinless binaries by tedious but straightforward calculations (the generalization of our Sec. 
HI through 3PN order is available in form of supplementary data; the more complicated 
case of spinning binaries is left for future research). Relying on the diploma thesis by one of 
the authors, we give explicit expressions for canonical center and relative coordinates 
in terms of standard canonical coordinates up to 2PN order. We particularly show that, 
using these expressions, the Pythagorean-theoreni-type Lorentz-invariant relation between 
the system’s total energy (or Hamiltonian squared), the rest energy (or mass squared, only 
depending on relative coordinates), and the total linear momentum squared holds up to 
2PN approximation. Working within a canonical framework we are always allowed to make 
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canonical transformations. This means that neither the Hamiltonian at any PN order nor 
the varions PN coordinate expressions we are calculating are uniquely dehned in the quasi- 
Minkowskian coordinates that we are using throughout our paper (recall the Schwarzschild 
metric in isotropic coordinates or in Schwarzschild coordinates but with cartesian compo¬ 
nents). Thus, no local physical interpretation of our PN expessions can be given. Readers 
interested in historical information on the relativistic two-body problem in relation to center 


and relative coordinates may hnd details in e.g. 


16]. 


II. FUNDAMENTALS 


We start with standard canonical position and linear momentum coordinates (variables) 
for the components of a binary system, denoted by (xa,pa), with coordinate-components 
(xl = Xai,Pai = pI) ^ usiug four-dimensioual metric signature +2, we do not need to dis¬ 
criminate between covariant and contravariant components of our cartesian coordinates -, 
where the labels a = (1, 2) and i = (1, 2, 3) number the involved objects and the coordinate- 
components, respectively. These variables are called the same in all our various PN ap¬ 
proximations. Also, the total linear momentum, P = total orbital angular 

momentum, L = ^ given by the same expressions through all PN orders. A 

simple transformation takes us to the following set of coordinates, which we call Newtonian 
coordinates, because in the realm of Newtonian physics they are nothing but the Newtonian 
center-of-niass, X^, the total linear momentum, P^(= P), the relative position, x^, and 
the relative momentum, p^, coordinates of two objects with mass parameters nia, 


X 

P 


N 

N 


mi 1712 

-^^-X2, 

mi -|- m2 mi -|- m 2 

Pi + P2, 


xN = Xi 


p^ = 


-X2, 

m 2 


mi -|- m2 


Pi 


mi 


mi -|- m2 


P2. 


( 1 ) 


Within the Hamiltonian framework a set of coordinates is called canonical if the equations 
of motion are in Hamiltonian form. 


X 


i 

a 


dH 

dpt 






dH 

dxt 


{PIH}, 


( 2 ) 


where {, } denotes the Poisson brackets. 

Every set of canonical coordinates can be used to evaluate the Poisson brackets and the 
canonical coordinates satisfy the canonical commutation relations, which read for the two 
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different sets of canonical coordinates, 


{Xai, Pbj] = Sab^ij , zero otherwise, 
{Xf, , [x^, p^] = 6ij , zero otherwise. 


(3) 

(4) 


We are looking for a new set of canonical coordinates (X, P, x, p) such that the equations 

of motion decouple into an external and an internal part. In Newtonian physics we simply 

2 2 

insert a Newtonian Hamiltonian {H = + I/(|xi — X 2 I)) and hnd that the equations 

decouple. 


* ^ J rn 


Pf = {Pf,P} = 


dH 


axN* 


T = 0 


x^,={xlH} = 

p^ = {pIh} = 


P^i 


dH 
dp^^ p 
dH 
dx^^ 


Gum ^ 

- 'T. 

' 12 


(5) 


where m = mi + m 2 , p = mim 2 /m, and ri 2 = |xi — X 2 I. Apparently, the quantities 
that determine the external motion (X'^, P^) are independent from those quantities of the 
internal motion (x'^, p^) and vice versa. 

In the relativistic theory we have to make an effort to decouple the equations of motion. 
We note that the energy-momentum relation, with c denoting the speed of light. 


P2 _ c2p2 = 


( 6 ) 


provides us with an equation, linking M, the invariant rest mass, and P, the conserved total 
energy. The total energy equals the Hamiltonian of the system since we are using canonical 
coordinates from the very beginning. It is a requirement for the notion of center and total 
momentum coordinates that M is not depending on external coordinates 


{W,M} =0, 


applying the canonical commutation relations 


{P,,M} =0, 


(7) 


{Xj, Pj} = 6ij , {xj, pj} = Sij , zero otherwise. 


( 8 ) 


The Hamiltonian determines the equation of motion of any function / on phase space, 

/ =lf,H} = + ( 9 ) 
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in particular, 


pk 



0 , 



( 10 ) 


( 11 ) 


Usually, in the binary problem, it is sufficient to just set = 0 in order to transform 
to the rest frame, e.g. In this paper though, we are interested in finding coordinates 
without resorting to this shortcut. We say that M is decoupling, if it is only depending on 
internal coordinates. 

Let us summarize the properties the new coordinates should have. First, the new set 
of coordinates has to be a canonical set of coordinates, so that the Hamiltonian form of 
the equations of motion is preserved. Second, the total linear and angular momenta have 
to be preserved in their forms and are therefore fixpoints of the transformation. Third, 
the equations of motion should decouple as described above. Last, the Newtonian limit 
should reproduce the known Newtonian two-body Hamiltonian in center-of-mass and relative 
coordinates, which is ensured because we do not leave the post-Newtonian framework. 

The first requirement is realised by introducing a canonical transformation Tg, which will 
raise the Newtonian coordinates to a given higher order post-Newtonian level: 


X = TgX.^ , X = T,x^, 

P = TgP^ , P = TgP^. (12) 


The transformation Tg can be represented in terms of its infinitesimal generator by a 
Lie-series of the form 

TgA = = A + {A, 9} + 2 9} , 9 } + q {{{^) 9 } , 9 } , 9 } + ■■ ■ ■ (13) 

Clearly, Tg written in this way, defines a canonical transformation. 

The second requirement then translates into the following conditions: 


{P, 9 } = 0, {J, 9 } = 0, (14) 

where J denotes the total angular momentum, i.e. J = L -|- S with S = [actually in 

our paper, the spins of the single objects are put equal to zero]. The generator g has 
been calculated in for various PN approximations. For the convenience of the reader, 

we recap the generators relevant to us in the Appendix. 
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The Poincare or inhomogeneous Lorentz group is the fundamental symmetry group 
of physical systems living in Minkowski spacetime. Isolated systems living in space- 
asymptotically flat spacetime, i.e. asymptotic Minkowski spacetime, also have the Poincare 
group as their symmetry group. This is the case here. However, our systems are not fully 
relativistic, but only relativistic to some given PN order. Therefore the Poincare group is 
not the exact symmetry group but only to a given PN order. In particular, this means that 
the Poincare algebra is only fulfilled to a given PN order as well. 

The Poincare algebra is explicitly given by (see ^ for details) 


{P\H} = {J\H}=0, 

{G‘, H} = P‘, {P‘,P>} = 0, 

{J‘,P’} = {/, Ji} = 

{J‘, G‘} = {G‘, P'} = 

{G‘, G‘} = 


(16) 

(16) 

(17) 

(18) 
(19) 


with P* = Pj, J* = = Gi, and where denote the total Hamiltonian, 

total linear momentum, total angular momentum including spin, and the center-of-energy 
position vector, respectively. For completeness, we like to note that the conserved Lorentz 
boost vector iP* is given by iP* = G* — fP*, where t denotes the asymptotic physical time 
(in flat spacetime, t would be a physical time parameter all-over space belonging to some 
specific inertial reference system). 


The relativistic canonical center position vector can then be written in the form, e.g. 


1^, 


X = 


Gc2 

H 


1 


M{H + Mc^) 


J - 


Gc2 

H 


X P X P 


This choice of X ensures that M will only depend on internal variables, i.e. 


DM 


- = {X\M} = 0, 


DM 


= {P\M] = 0. 


( 20 ) 


( 21 ) 


aP* ' ’ ' ’ dX^ 

The Poincare algebra also indicates the following Poisson brackets, which we want to add 
for completeness. 


{M, P} = 0, ^{X, P} = P. 

In our case of non-spinning objects, the canonical center position vector simplifies to 


X = 


Gc2 

H 


M {H + Mc2) 


L- 


Gc^ 

H 


X P X P. 


( 22 ) 


(23) 
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Notice that the canonical center position vector does not coincide in general with the center- 


of-energy (or center-of-mass) position vector which is defined by e.g. 
of the varions center definitions see e.g. H- 


17| . For a synopsis 


III. CANONICAL CENTER AND RELATIVE COORDINATES 


With the aid of the relations in equation (1121) and the generator g derivec 


can now ca. 
XTensor 


culate the new coordinates explicitly. We use the Mathematica 


20 


in 


141 ■ we 


19j package 


21| . All files are available upon request. 


The coordinate transformations are expressions with many summands, therefore, we 
split them into terms proportional to First, we start with the center coordinate X, 

which, initially, is identical to the center-of-mass X^. Second, we present the new relative 
coordinate x and last, the new relative momentum p. The total linear momentum is left 
unchanged, i.e. P = = Pi + P 2 - 

The new coordinates are functions of (P^,x'^,p^) or, after a simple substitution, 
(x'^,Pi,P 2 ). Apart from X, there is no dependency from X^ by construction. We use 
the abbreviation rj 2 = 


results for X 


X* = X"" + + XjpN 4- 0(c-“) 


(24) 


-yk _ 

^IPN — 2 


{x'l - 4) 


^ {mi - m2)(pi ■ P 2 ) Gmim2{mi - m 2 ) 
2m2m^ 2mim^ 2m^ 2ri2m‘^ 


{m2p\ — mip\) 
2m^ 


(Pi + P 2 ) ■ (xi - X 2 ; 


-yk 

"^2PN 


p\ \ (pi • (xi -X 2 ))( - 


Gmi{mi — 5m2)m2 
4 (mi 7712)^12 
3m2(mi 2m2)(pi)^ 3m2(pi ■ P 2 ) 


8mi(mi-|-m2)^ 4(mi-|-m 2 )® 


-2m1 -|- 2mim2 + 7n|)(p2)^ 
8m2(mi -|- m 2 )® 
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+ (p2 ■ (Xi - X2)) 


G{5ml + 16m^m2 + + Sm^) 

4(mi + m2)Vi2 


3 m 2 (mi + 2 m 2 )(pi)^ 3 m 2 (pi ■ P2) (- 2 m^ + 2 mim 2 + m|)(p2)" 


+ 


+ P2 


(Pl ■ (xi - X2 


8mi(mi+m2)^ 4(mi + 7712 )^ 8m2(mi+m2)® 

G(5mf + 10m^m2 + + Sm^) 


4(mi + m2)^ri2 


{mf + 2mim2 — 27712) (pi)^ 


87771 ( 7771 + 7772 )® 4(7771 + 7772 )® 

/ G 777 i 7772 (-57771 + 7772) 

(777^ + 277717772 - 27773) (pi)^ 3777i(pi-p2) 


37771 (pi • P 2 ) 37771 (27771 + 7772) (p2)^ 


87771(7771 + 7772 )® 


+ 


87772(7771 + 7772 )® 


37771 (27771 + " 2 - 2 ) (P 2 )' 


4 ( 7771 + 7772 )® 87772 ( 7771 + 7772 )® 


(4 


-4) 


G^777i(777i — 7772)7772(777^ + TTl^) 7773 ( 777 ^ + 5777i7772 + TTl^) 


4(7771 + ?T72)®?"i' 


3(7771 — 7772 ) 


I 2 8777 ^( 7771 + 7772 )® 

Gmi{ 5 ml + Qm\m2 + Qmim^ — m\ 


(Pi)^ 


4(7771 + ^^ 2 )® 
m\{m\ + 5777 i 7772 + 7773) 
87773(7771 + 7772)® 

( 777 ^ + 5777i7772 — 87773 ) 


/ \2 Vjr// \JIIL\IIL<2 ^2f/ / \2 

(Pi ■ P2) H- Tim. —r+r+z:-(P2) 


47772(7771 + m2)^rx2 
4 Gm 2 { jn \ — 677747772 — 677717772 — StTI^) 2 
47771(7771 + 7772 )^ri 2 


87771(7771 + 7772)® 


(P.) = (P. ■ P 2 ) + 5 ”" , "''I (Pl)"(P 2 ) = 


G(mi — m 2 ) {7ml + 877717772 + 7ml) 

4(7771 + ■'272)^12 
(— 8777 ^ + 5777i 7772 + TTl^) 


8(7771 + 7772 )® 

(Pl ■ P2) 


87772(7771 + 7772 )® 

G(777i — m 2 ) 


(P2) (Pl • P2) 


4(7771 + m2)'^r\2 


Pl ■ (xi - X2))(P2 ■ (xi - X2)) 


(25) 


B. results for x 


= X^k X 4 PN + XgpN + 0{c ®) 


(26) 


fc _ I f (-7771 + 27772 ) k {‘2mi - m 2 ) 
^IPN ^2 ^27771 (7771 + 7772)^^^ 27772(7771 + m2) 


rP2 


(Pl + P2) ■ X 


,N 


(27) 





























^2PN — 




P 2 ' 


„Nfc 


X 


G(7mf — 4m^m2 — 9mim| — 10m|) 

4mi(mi + m2)^ri2 

{—3ml + 4mfm2 + Smim^ + 4m|)(pi)^ 


+ 


8mf (mi + m2)^ 

(—2mf — m\m 2 + 2mim2 + 4m2)(pi • P 2 ) 
4mim2(mi + m2)^ 

(4mi + 3mim2 — 4m|)(p2)^ 


(pi + P 2 ) ■ X 


N 


+ 


8mim2(mi + m2)^ 

G(10mf + 9m\m2 + 4mim2 — Tm^) 
4m2(mi + m2)^ri2 
{—Ami + 3mim2 + 4m2)(pi)^ 
8mim2(mi + m2)^ 

(4mi + 2mlm2 — mim^ — 2ml){pi ■ P 2 ) 

4mim2(mi + m2)"^ 

(4mi + 8mlm2 + 4mim| — 3ml){p2)‘^ 


5G 


8{mi + m2)ri2 
G{mi — m 2 ) 


8 {mi + m2)^r’i2 


8m2(mi + m2)'^ 

(Pl +P2)^ 

(mi - 5m2)(pi ■ x’^) 


,N 


^Pl + P2j ■ X 


+ (5mi - m 2 )(p 2 ■ x’^) j (^(pi + P 2 ) ■ 


(28) 


C. results for p 


P*=P"*+P?PN+PiN + C’(c-“) 


(29) 


PlPN 


/ Gmim2(-mi + m 2 ) 


(P? + P 2 ) 


m|(pi 


2(mi + m 2 )^ri 2 


^ (-mi + m 2 )(pi ■P 2 ) ^ 
2(mi + m2)^ 


2mi(mi + m2)^ 

mKP2)^ 

2m2(mi + m2)^ 


+ " 2(mi + «P‘ + K) ■ " ) 


(30) 
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P2PN 


(K + 


I — m2)m2{rn\ + m^) m|(mi + 5mini2 + ?Ti2)(pi)'^ 


+ 


4(mi + m2)^T’i2 8mf(mi + 7712 )^ 

3(mi - m2)(pi ■ P2)^ ml{ml + 5mim2 + m|)(p2)'^ 


+ (Pi 


4(mi + m 2 )® 8m|(mi+m2)® 

^2(ijn\ + 5mim2 - 8m2)(pi ■ P2) 3 (mi - m2)(p2)^ 
8mi(mi+m2)® 8(mi+m2)® 

(—8m^ + 5mim2 + m2)(p2)^' 




(Pi P2) ^ 8m2(mi+m2)® 

G(10mf + 12 mfm 2 + ITm^m^ + 8 mim 2 + 5 m 2 )(p 2 )^ 
8m2(mi + m2)^ri2 

Gm 2 (—7mf + 2m\m2 + 7mim\ + lOm®) ^2 

■(Pi) 


8mi(mi + m2)Vi2 

G(7mf — 4 m^m 2 — llmim^ — 12m|) 

4(mi + m2)^ri2 

Gm2(-mi + m2)(pi ■ G(mi - m2)(mi + 2m2)(pi • x®^)(p2 ■ x®') 


(Pi • P2) 


+ 


pV 


8(mi + m2)®r^2 

3G(mi - m2)m2(p2 • x^)^ 

8 (mi + m2)®r’i2 

Gmi(10m5 + 7mlm2 + 2mim2 — 7m2)(p2)^ 

8m2(mi + m2)'^r'i2 

G(5m^ + 8 mfm 2 + 17mim2 + 12mim| + lOm^) 

8mi(mi + m2)Vi2 

G(12mi + llmim2 + 4mim| — 7m2) 

4(mi + m 2 ) V 12 

,N\2 


4(mi + m 2 )®ri 2 


(Pi)^ 


(Pi ■ P2) 


3Gmi(mi — m2)(pi ■ x®^)^ G(mi — m2)(2mi + m2)(pi ■ x®^)(p2 • x^) 


8 (mi + m2)®r’i2 
Gmi(-mi + m2)(p2 ■ x^)^ 


+ 


4(mi + m2)®ri2 


8 (mi + m 2 )®rf. 


12 


+ X 


Nk 


3G(mi — m 2 )^m 2 (pi ■ x^)® G^mi(mi — m 2 )m 2 {jn\ + m 2 )(p 2 


.N'l 


8(mi + m2)'^7’i 


12 2(mi + m2)®rf2 

3G(mi — m 2 )( 3 mf + mim 2 + 4m|)(pi ■ x'^)^(p 2 • x^) 
8 (mi + m 2 )^ri 2 
3Gmi(mi — m2)^(p2 ■ x^)® 


8 (mi + m 2 )^ri 2 


(Pi ■ P2) 


5G(mi — 2m2)(pi • x®^) 5G(—2mi + m 2 )(p 2 • x^) 


,N^ 


4(mi + m2)^pf. 


+ 


12 


4(mi + m2Yr\. 
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, + 6m^m2 + 5mim2 + 5nil){pi ■ x^) 

(P 2 ) - 


8m2(mi + m2)^r^2 
Gmi{10m‘l + mim2 — 5m2)(p2 ■ x^) 

8m2(mi + m2)^rf2 

Gm2{—5ml + mim2 + 10m2)(pi ■ x^) 

8mi(mi + m2Yr\2 

G{5ml + bm\m2 + Qmim^ + 10m2)(p2 • x^^ 


(pir - 


(pi -x^)! - 


8mi(mi + m2Yr\2 
G‘^mi{mi — m2)m2{m\ + m|) 


2(mi + m2)^r^2 
3G(mi — m2)(4m^ + mim2 + 3m2)(p2 ■ x^)^ 
8(mi + 7712)^^2 


(31) 


The total angular momentum of our binary system with non-spinning components, i.e. 
J = L = J2a^a X Pa, in Center and relative coordinates, reads. 


X + x^ X p^ = X X P + X X p = L. 


(32) 


We confirmed that J = L is a fixpoint of our canonical transformation. 


IV. INVERSE TRANSFORMATION THROUGH 2PN 

It is necessary to know the original Newtonian canonical coordinates in terms of the new 
center and relative coordinates, in order to substitute the old for the new coordinates in 
a given equation. With a simple further transformation, the standard coordinates for the 
single components of the binary system can then be obtained as functions of the new center 
and relative coordinates. The inverse transformation through 2PN order reads, 

x’^ = X - XIPN - X 2 PN + G(c“®), (33a) 

p’^ = p - PiPN - P2PN + c>(c"®), (33b) 

X^ = X-XiPN-X 2 PN + C»(c-®). (33c) 

The equations (l33|) are simply rearrangements of previous relations and the single terms 
xiPN, X 2 PN, PiPN, P 2 PN, XiPN, X 2 PN are still functions of (x^,pi,p 2 ). Now we are going to 
invert the transformations step by step. In a first step we replace simultaneously x^ with 
X and p'^ with p on the right hand side of equations fl55]) . Then we regroup the terms 
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proportional to c ^ and find the inverse transformation to first order. The hrst order inverse 
transformation is then used in another step to derive the second order inverse transformation. 




pk 


fP-x) 




pk 


2(mi + 7712 )^ 

_^_( 

8(mi + m2Y 

1 


P 


- m 2 )(P ■ x) 
m\m2 + rriim^ 


P)"(P-x) + 


(p)^(P-x) + G 


(-mi + m 2 ) 
2mim2(mi + m2Y 

(5mf + 6mim2 + Sm^ 


P 


+ 


X 


2mim2(mi+m2)^'^' ' A{ml+m2Y^/l 

f. I (—mi + m2)(P)^(P • x) (mf — mim2 + m^) 


P P)(P x) 


fP-x) 


X ■ X 


2mim2(mi + 1 x 12 )^ 


+ 


mim2(mi + m2)^ 


P-p)(P-x) 


(—m^ + ml) \ I ry (~5mf — 3m|m2 + 3mim| + 5m|) 


2m\m\ 


-(p)^(P-x) + G 


-G 


_ 3(mi -m2)^(P -x)^ 
8(mi + m2)^(x ■ x)3/2 

5(P)2 


2mim2(mi + m2)^-\/(x ■ x) 

'-mi + m2)(P • x)(p ■ x) 
2(mi + m2)^(x • x)3/2 


P x) 


8(mi + m2) v^(x ■ X) 


(34) 
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+ 


p 


(P-P) 


+ 


(—mi + m2)(p)^ Grriiimi — m2)m2 


+ 


+ X 


2(mi+7712)2 2mi7?72(mi + 7772) 2(7?7i + 7772)2 a/C 

f. Gmim 2 {—mi + 7772)(P ■ x 


X • X 


2(7771 + 7T72)2(x ■ x)^/^ 


pk 


(P)^(P-P) 

8(7771 + 7772 )^ 


+ 


(7771 - 7772) (7771 + 7772) (p)^ 


(P)^(P-P) 


8777i7772 


G(5777i + 6777i7772 + 57772 
4(7771 + 7772 )^-^/(x • x) 
G(777 i — 7772)2(P ■ x)(p ■ x) 


(P ■ P) + 


277717772(7771 + 7772)2 

G777 i 7772(—777i + 7772)(P ' x)2 


+ 


4(7771 + 7772)'^(x ■ x)3/2 
G(777i - 7772 ) (p ■ X)2 


4(7771 + 7772)^(x ■ x)3/2 4(7771 + 7?72)2 (x ■ x)^/^ 

(727771(7771 — 7772)7772(7772 + 7772) ^ G{mi — m 2 ){ 5 m‘l + 877717772 + 57772 )(p )2 


4(7771 + 7?72)^(x ■ x) 


477717772 (7771 + m2y ^/l 


X • X 


p 


X 


5c;(p)^ 


G{mi — 7772)2(P ■ x)2 G{mi — 7772)(P ' x)(p ■ x) 

8(7771 + 7772)3(x • x) 3/2 ^ 2(777i + 7772)2(x • x) 3/2 ^2) 

f 3(7(7771 — 7772)^(P ' x) 2 (p • x) 3(7(—777i + 7772)(P • x)(p ■ X 


X ■ X 


I 8(7771 + 7772 )^(x • x)®/2 

(727771(7771 — 7772)7772(7772 + 7773) 


+ 


4(7771 + 7772)2(x ■ x)^/2 

G{3mi + m 2 ) {mi + 87772 ) 


P-X + 


2(7771 + 7772 )^(x ■ x)2 
G{mi — 7772)(5777i + 8777i7772 + StTI^) 
477717772(7771 + 7772)2 (x ■ x)3/2 

Gmi{mi - 7772)7772 ^p ^2 _ 


4(7771 + 7772)^(x • x )^/2 


(P-p)(P-x) 


(p)2(P.x) 


5C; 


4(7771 + 7?72)^(x • x)3/2 


8(7771 + 7772 ) (x • x)3/2 


:p-x)(p) 


(35) 
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1 

ri2 


VI 


X ■ X 


(P ■ x )2 ^ {-mi + m2)(P ■ x)(p ■ x) 

2{mi + m2)^(x • x)^/2 mim2{mi + m2)(x ■ x)^/^ 

3 (P ■ x)^ 3 (mi — m2)(P ■ x)^(p • x) 


3 (mi — m2)^(P ■ x)^(p • x)^ 


8(mi + m2)‘^(x • x)^/2 2mim2(mi + m2)^(x ■ x)®/^ 2m^m2(?7ii + m2)^(x • x)^/^ 

G{7m\ + 18 mim 2 + 7m2)(P ■ x)^ ^ G{mi — m2)(5mf + 9mim2 + 5m2)(P ■ x)(p • x) 


8(mi + m2)^(x • x)^ 


2mim2{mi + 7712)^1 


X ■ X 


(P)2(P.X) 


(mi - m2)(P • p)(P ■ x) 


, , , ^ - _ (mf + m|)(p)^(P • 

2(mi + m2)^(x • x)3/2 mim2(mi + m2)^(x ■ x)^/^ 2mfm2(mi + m2)^(x 
(mi - m2)(P)^(P ■ x)(p ■ x) _ (mi - mim2 + ml){P ■ p) 

2mim2(mi + m2)^(x ■ x)^/2 rnfm^mi + m2)^(x ■ x)^/^ ^ 


(mi - m 2 )(mi + m2)(p)^ 
2mim2(x ■ x)3/2 


;p-x)(p-x) + 


5 G'(P)^ 


8(mi + m2)(x ■ x) 


• x)3/2 


( 36 ) 




fp-x) 


-p 


2(mi + m2)2 

(P-P) 


(mi—m2)(p)^ Gmim2(—mi + m 2 ) 

~t~ o ^ 77 o 


2 (mi + m2)2 2mfm2 + 2mim| 2(mi+m2)2v^ 


X 


X ■ X 


P 


3 (P)^ 


^ (—mi + m2)(P • p) ^ (mi — mim2 + m|)(p)^ 


8(mi + m2)^ mim2(mi+m2)3 
G( 3 mi + m2)(mi + 3m2)(P • x 


2mfm2(mi + m2)^ 


4 (mi + m 2 )^Vl 


X ■ X 


F' 


(P-p) 


+ 


(-mi + m2)(p)' 


(2(mi + m2)'^ ' 4mim2(mi + m2)3 

Gmi(mi — m2)m2(P • x) 5 G(p • x 

4 (mi + m2Y \/(x ■ x) 4 (mi + m2) VI 

p)2 


P x) 


IX-XJ 

,2 I 


X ■ X 


X 


{ -mj + ml) {pY 
8m3m| 


(p p: 


+ 


ip) 


8(mi + m2)^ 2mim2(mi + m2)^ 


fP-x) 


Gmi(mi — m2)m2(P • x)^ G(mi — m2)^(P • x)(p • x) G(—mi + m2)(p ■ x)^ 
2(mi + m2)^(x ■ x)3/2 4(mi + m2)3(x • x)3/2 4(mi + m2)^(x • x)3/2 

G^mi(mi — m 2 )m 2 {m\ + m^) ^ G( 5 mi + 6mim2 + 5 m 2 )(P ■ p) 


4 (mi + m2)^(x ■ x) 4(mi+m2)^yC 

G(mi — m2)(5mi + 8mim2 + 5m2)(p)^ 


X ■ X 


4mim2(mi + m2YVi 


X ■ X 


(37) 
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We thus found the old coordinates as functions of the new coordinates. Notice the 


tremendous simplification of the relations in the rest frame where P = 0 holds. 


V. DECOUPLING THE REST MASS 


In this last section, we want to show that our newly derived coordinates do as promised, 
i.e. M decouples in the new coordinates up to 2PN approximation, analogously to the 
Newtonian case. For convenience, let us restate equation ([6]); 


r .4 ^2 


(38) 


After squaring the Hamiltonian H = mc^ + + ^iPipN + ^hf 2 PN + 0{c~^) and grouping 

all terms together that are proportional to a factor in we find the following structure: 


= Bo + + 4^2 + 4^3 + 0(c-*). 


(39) 


In the spirit of the perturbative nature of post-Newtonian theory, we expand M with the 
parameter in a Taylor series: 


c rp2 C 
-^0 


and find, 

M = mi + m2 + 


1 El 1 AE 2 EQ — (Ei)'^ 1 Ef — AE 1 E 2 EQ SE^Eq 


+ 


8E^ 


16E^ 


+ 0{c 


-83 


(40) 


Gmim2 ^ {mi + m 2 ) 


ri2 


+ 


G‘^mim2 {mi + m 2 ) 


2mim2 

1 


(41) 


2rf2 


Um\ml{V^ ■ 


Smfml {mi + m2) \ 

Amim2 {—mi^ + rn2^) (P^ ■ p^)(p'^)^ + {m\ + m\m2 + mim\ + m^) (p^)' 
G(mim2 (ml — ml') ri2(P’^ ■ p'^) 


2mim2 {mi + m 2 ) 

-I- (rui -I- 7712 )^ (3ml + 7mim2 + 3ml) r^2(p'^)^ 
-|- mim2 ^ — mim2{P^ ■ 


+ {ml - ml){P^ ■ x^)(pN ■ x^) + {mi + m2)\p^ ■ x^)^)) 


It is obvious that M does not decouple in the old Newtonian coordinates. 
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We now substitute for the new coordinates in M: 


M 



1 (^1 + ml) 4 G' 

Smfml 


2(x • x) 


G 





2mim2(x ■ x)^/2 
1 (mf + ml) , f 


g G^mim2('^^ ^'^-m2 + m|) 



^2 ("^1 + ^ 


2mim2(x 



G 


(^ - 3m^m^(p ■ x)^ - 2m?m^(p)^(p ■ x)2(x ■ x) 
13m^m2 + 5m2)(p)'^(x ■ x)^j . 


Brnfm^x ■ x)^/2 


(42) 


We see, very nicely, that M, in our newly derived coordinates, solely depends on relative 
variables and thus decouples. 
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Appendix A: PN generators 


Rothe and Schafer derived the generator g for the obtention of the canonical center and 
relative coordinates in denoted therein (?point (“point” meaning without spin). Also, the 
abbreviation ni 2 = (xi — X 2 )/ri 2 is used. 


fi'point — - [S'lPN + fi'2PN + (1 Ut 2)] 


(Al) 


giPN — c 


-2 


rrin 


2m' 


mim- 
— m 


(ni2 ■ Pi)Piri2 + 


ms 


mim-^’ 




(ni2 ■ pi) (pi ■ P 2 ) ri2 


■ (ni2 ■ P 2 ) pjru 

/i( 2 m 2 - m) 

G -(ni2 ■ Pi 

m 


(A2a) 
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5'2PN 


^-4 


ml{m2 + 2mi) . x / 2 n 2 m|(m2 + 2mi) ^ / 2\2 

4m^m‘ ' P‘> - 4mfm‘ 


2m2 + /i — m 


I - ' ( \ i ^ 2 I 2m2 + /i?77, 2 

H-^^5—5— (ni2 ■ Pi) (Pi ■ P2) Pir’12 H-^— (ni2 ■ P2) (Pi ■ P2) ^1^12 


2mf rr? 

2m2 — m , 

-^-( 

2/i 

\ml{m2 - mi) 

- :—~A - 

4m^ 

2 /i(3/x — 8m) — 5m2(/i — 2m) 

4m im^ 

6/i^ — 5m2/x — 16m/x + 5m^ + 5 m 2 m 


2mf m^ 

2 2m2 — m \ 2 2 

ni2 ■ Pi) (Pi ■ P2) ri2---— (ni2 ■ Pi)P2Pin2 

4/i m^ 

3 mi(-5m2 + 6/i + m) 2 

(ni2 ■ pi) H-—^-(ni2 ■ P2) (ni2 ■ Pi) 


4m im^ 

6m^ — 17 m2m — 3yum + 6m2/U 


4m^ 
ni2 ■ Pi) Pi 

(ni2 ■ P 2 ) Pi 


2m^ 


(ni2 ■ pi) (pi ■ P2) 


+G' 


2 p( 2 yU - m)(m - 2 m 2 ) (ni 2 ■ pi) 


2 m 


?^12 


(A2b) 
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